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Abstract 



Presented here is a study of the large pt J/ ifr production in hard diffractive 
process by the pomeron exchange at the Fermilab Tevatron. We find that 



^ ■ this process (pp — > p + J/ip + -X") can be used to probe the gluon content 



of the pomeron and to measure the gluon fraction of the pomeron. And the 
diffractive direct J/ip production can also provide another crucial test for the 
color-octet fragmentation mechanism. Using the renormalised pomeron flux 
factor D ~ 1/9, the single diffractive J ftp production cross section at large 
Pt (> 8GeV) is found to be of order of O.Olnb, and the ratio of the single 
diffractive to the non-diffractive J/ip production is 0.65 ±0.15% for the gluon 
fraction f g = 0.7 ± 0.2. 
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In the past few years there has been a renaissance of interest in diffractive scattering. 
Diffractive processes in hadron collisions are well described by the Regge theory in terms of 
the pomeron (IP) exchange The pomeron carries quantum numbers of the vacuum, 

so it is a colorless entity in QCD language, which may lead to the "rapidity gap" events 
in experiments However, the nature of pomeron and its reaction with hadrons remain a 
mystery. In |§, the hard diffractive scattering processes have been suggested to resolve the 
quark and gluon content in the pomeron. That is to say, the pomeron has partonic structure, 
just as hadrons and nuclei. Therefore, various processes may be considered to probe the 
partonic structure of the pomeron at high energy colliders . 

On the experimental side, the UA8 collaboration at CERN SppS collider have studied 
diffractive dijet production at yfs = 630GeV [[/], which indicate a dominant hard partonic 
structure of the pomeron. The HI and ZEUS have studied the diffractive deep inelastic 
scattering (DDIS) and dijet photoproduction in ep collision at y/s = 300GeV ||. From 
these measurements, ZEUS determined that the gluon fraction of the pomeron f g is in the 
range of 0.3 < f g < 0.8, and HI determined the quark fraction of the pomeron f q « 0.2. The 
partonic structure of the pomeron is also studied recently by the CDF collaboration at the 
Fermilab Tevatron through the diffractive W production || and dijet production |TD|, which 



give further evidence for the hard partonic structure of the pomeron. The combination of 
these two measurements determined the gluon fraction of the pomeron to be 0.7 ± 0.4. 

In this paper, we will discuss another diffractive process, the single diffractive (SD) J/ip 
production at large pt, (shown in Fig.l), 

p + p^p+J/ip + X. (1) 

By the following calculations, we will show that the SD J/ip production is sensitive to 
the gluon fraction of the pomeron. So, the measurement of diffractive J/ip production at 
the Fermilab Tevatron would provide a probe of the gluon distribution in the pomeron. 
Importantly, at hadron colliders the J/ip production is of special significance because it 
has extremely clean signature through its leptonic decay modes. Furthermore, the SD J/ip 



production is also interesting to the study of heavy quarkonium production mechanism, 
which is another hot topic in the past few years. 

For a long time, it was believed that the heavy quarkonium production at large px 
dominantly comes from the leading order color-singlet processes ||11|| . But, as pointed out 
by Braaten and Yuan [12], the fragmentation contributions may dominate over those from 



leading-order processes at sufficiently large pt, although the fragmentation processes are of 
higher order in strong coupling constant a s . However, the measurements of large pt char- 
monia production from the CDF at the Tevatron show a large excess of direct production 
(excluding the contribution from b decays and the feeddown from Xc) both for J/ip and ip' 
|T3| fT4|| . The experimental measurement is a factor of 30 ~ 50 larger than the theoreti- 
cal prediction of the Color-Singlet Model even if including the fragmentation contributions. 
Motivated by this "surplus" problem, a new mechanism for heavy quarkonium production at 
large pt in hadronic collisions, named as Color-Octet gluon fragmentation has been proposed 



I5| , which is based on the factorisation formalism of non-relativistic quantum chromody- 
namics (NRQCD) ]nj. In the past few years, applications of the NRQCD factorisation 
formalism to J/i/j(i[)') production at various experimental facilities have been studied (JTT . 

According to NRQCD factorisation formalism, the gluon fragmentation to J /if) produc- 
tion can be factorized as 

D^j/^z^ 2 ) = Y,d g ^ n (z,» 2 ){0 J n ^), (2) 

n 

where z is the longitudinal momentum fraction carried by the produced J/ip in gluon frag- 
mentation, /i = 2m c is the fragmentation scale. d g ^ n represent the short-distance coefficients 
associated with the perturbative subprocesses in which a cc pair is produced in a configu- 
ration denoted by n (angular momentum 2S+1 Lj and color index 1 or 8). (O^^) are the 
long distance nonperturbative matrix elements demonstrating the probability of a cc pair 
evolving into the physical state J/ijj. The short-distance coefficients d g ^ n can be obtained 
from perturbative calculations in powers of coupling constant a s . (O^J^) consist of two 
kinds of matrix elements, i.e., the color-singlet and color-octet matrix elements (according 



to the color index is 1 or 8). 

For J/ip production in gluon fragmentation, the color-octet matrix element (Og^( 3 Si)} 
is smaller than the color-singlet matrix element (0(^( 3 Si)} by a factor of order v 4 according 
to the NRQCD velocity scaling rules. However, the short-distance coefficient for the color- 
octet term in Eq.(0) is larger than that for the color-singlet term by a factor of order 
l/ojg. Numerical results show that color-octet contributions are 50 times larger than color- 
singlet contributions |15|| . In the following calculations, we neglect the color-singlet term in 
gluon fragmentation in Eq.@, and only consider the color-octet gluon fragmentation. The 



leading-order color-octet gluon fragmentation to J / ip production gives |L5 

In our calculations, the effects of the evolution of gluon fragmentation function with scale 
fi 2 are neglected, which may introduce some error. However, as argued in [18], including 



evolution would not necessarily be an improvement, since naive Altarelli-Parisi equations do 
not respect the phase-space constraint D g ^j/^(z, n 2 ) = for z < m 2 ,^//! 2 [[TJJ. 

As shown in Fig.l and Fig.2, the SD process pp — > p + J /if} + X by pomeron exchange 
consists of three steps. First, a pomeron is emitted from the proton with a small squared 
momentum transfer, t = (pi — Pf) 2 , where pi and pf are the momenta of the initial and the 
final states of the proton respectively. Second, partons interaction between the pomeron 
and the antiproton takes place in the large momentum transfer processes (see Fig.2). In 
the third step, J/if> is produced via the fragmentation processes. Because the fragmentation 
contributions dominantly come from the color-octet gluon fragmentation process, in our 
calculations we calculate the gg(qq) — > gg, q(q)g — > q(q)g processes and followed by the 
gluon fragmentation Eq.(Q). 

Using the pomeron factorisation formalism [0,f§], we write the SD J /if) production cross 
section as 

da(pp — > p + J /if) + X) 



d^dt 



fjP /p (U)dc7(lPp^J/iP + X), (4) 



where £ is the momentum fraction of the proton carried by the pomeron. fjp/ p is the pomeron 
"flux" factor, 



jp/p\s,t) = — jpp, , - = -r^r? ^ \t) 



a£ p {s>,t) 16tt 
K^- 2a ^F 2 {t). (5) 



Following [22 1, the parameters are choosen as 



K = 0.73GeV~ 2 , a(t) = 1 + 0.115 + 0.26{GeV~ 2 )t, F 2 (t)=e 4M . (6) 



In our calculations, we use the renormalizedfiux factor £)•/(£, t) [22], which may preserve 
the shapes of M 2 and t distributions in SD and predict the experimentally observed SD cross 
section at all energies. Here D is defined as 

D = Mm(l, l/N), (7) 

with 

£,max OO 

M= J di J fr/ P (U), (8) 

(mm ^ 

where, £ m j n = M^/s with Mq = 1.5GeV 2 (effective threshold) and £ max = 0.1 (coherence 
limit). For the SD process at the Tevatron (-y/s = lSOOGeV), the renormalized factor 
D w 1/9. As a conservative estimate, we will take this as a tentative value for D in the 
following calculations. (However, we should keep in mind that this flux factor D has not 
been well determined experimentally. If the precise value can be obtained in the future, our 
results will change accordingly.) 

On the parton structure functions of the pomeron, we assume the hard form [|B|P,|T0| , 



PG((3) = (f q + f g ) [6/3(1 -0)], (9) 

where f3 is the momentum fraction of the pomeron carried by the quarks and gluons. f q and 
f g are the quark and gluon fractions of the pomeron respectively. The momentum sum rule 
constrains f q + f g = 1. We neglect any Q 2 evolution in the above parton densities of the 
pomeron [|J. 



With these partons densities in the pomeron and the parton distribution functions in 
the antiproton, daiJPp —>■ J/ip + X) can be calculated by employing the usual way in the 
parton model calculations in hadronic collisions. We use the MRS (A) parton distribution 



functions |2T| to generate the production cross section, and set the renomalization scale and 
the factorization scale both equal to the transverse momentum of the fragmenting gluon 
A 4 = Pr{g) ~ Pt(V')- m gluon fragmentation, the input parameters are taken to be 

m c = 1.5GeV, a s (2m c ) = 0.26, (Og 7 ^ 3 ^)) = OMOQGeV 3 . (10) 

The value of the color-octet matrix element (Og^ ( 3 Si)) follows the fitted value in |2(J by 



comparing the theoretical prediction to the experimental data at the Tevatron. In the total 
cross section, we also include the contributions from the Xc and ip' feeddowns through gluon 
fragmentation g — > Xc and g — > ip' followed by Xc - > ^7 an d ip' ipX. The feeddown 
contributions give the same pr distribution of J/ip and contribute about one third of the 
total prompt J/tp production cross section (see Ref. flUf ). The leptonic decay branching ratio 
Br(J/ip — > = 0.0597 is also multiplied in the cross section. The above described 

procedure (including the fragmentation approximation) can reproduce the large px J/if > 
production in the central region (i.e, \rj(J/ip)\ < 0.6) at the Tevatron [ TA ,|T5|] . 



Because the qq annihilation process only contributes a small portion to the total cross 
section of large pt J/ip production, the SD J/tp production is insensitive to the quark 
flavor number of the pomeron. In our calculations we only consider two-quark flavors. A 
pseudorapidity cut of 2.0 < rj(ip) < 4.0 was also performed on the produced J/ip. The 
diffractive variable £ and t are integrated over the range of < £ < 0.1 and \t\ < lGeV 2 . 
In Fig. 3, We show the cross section of SD J/ip production as a function of the minimum 
transfer momentum of the produced J/ip. Because the gg process is the dominant process 
in the production of J/ip, the cross section is sensitive to the gluon component in the 
pomeron (i.e., the gluon fraction f g ). In this figure, we plot three curves correspondence 
to three different values of the gluon fraction f g . The solid line represents the cross section 
a(pp — > p + J/i/j + X) x Br(J/tp — * for the gluon fraction set to be f g = 0.7 (which is 



determined by the experiments at the Tevatron HlPp . The dotted line is for f g = (f g = 1) 
and the dashed line for f g = 1 (f q = 0). These curves show that the SD J/i/j production 
cross section may reach the level of order of O.Olnb (for pT(min) = 8GeV), and therefore is 
observable at the Fermilab Tevatron at present. 

In the above calculations, we use the widely used parametrization of the gluon distri- 
bution in the pomeron, i.e., the hard form Eq.(|9|). However, the precise form of the gluon 
density of the pomeron is unknown at present and this will affect the px distribution of the 
SD J/ijj production. Different parametrizations will give rise to different spectra. If the 
gluon in the pomeron is soft, e.g., the gluon distribution in the pomeron behaves like, 

/3G(/3)=6(1-/?) 5 , (11) 

the spectra would be different. In Fig.4, we show the px distributions of the cross section 
of the SD J/ip production in both the hard gluon and soft gluon cases. The result shown 
in Fig.4 is consistent with the expectation that the softer gluon will favor J/ip production 
with smaller px, while the harder gluon will favor larger px- The differential cross sections 
at large px for these two cases are different, but their differences are not so critical to 
distinguish between them. So, in the following discussion, we mainly limit us to the hard 
gluon parametrization of the pomeron, but we will also mention the result for the soft gluon 
parameterization. 

As a comparison, we also calculate the non-diffractive (ND) forward J/ip production 
(p + p — > J/ip + X) in the same kinematic region, i.e., 2.0 < f](ip) < 4.0. The forward region 
J/ip production is also interesting to the study of the J/tj) production mechanism, because 
the relative contributions of different mechanisms may vary with J/ijj rapidity. Furthermore, 
the comparison of forward and central region J j ip production can provide a consistent test 
of the "surplus" problem of ip' and J/ip found at the Tevatron ||14[| . In Fig. 5, we plot the ND 
forward J/ip production cross section as a function of px(min). This theoretical prediction 
can be used to compare with the experimental data, and may provide a further test for 
the color-octet gluon fragmentation production mechanism. (The DO collaboration at the 
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Fermilab Tevatron have reported the forward J / ip production data [ 24[] . However, these data 
do not exclude the contributions from b decays. We hope that the prompt J/ip production 
data in the forward region may be obtained in the near future.) 

In Fig. 6, we plot the ratio R{ip) = o'sd/o'nd as a function of f g (solid line), where asn is 
the cross section for SD J/ip production, and cttvd is for the ND J/ip production in the same 
kinematic region. The kinematic constraints for the produced J/ip are the same for these 
two processes, i.e., £>T(min) = 8GeV and 2.0 < t](ip) < 4.0. The ratio R(ip) increase from 
0.1% as f g = to 0.9% as f g = 1.0. For f g = 0.7 ±0.2, the ratio R(ip) will be 0.65 ±0.15%. 
We must note that R(ip) is independent of the choice of the color-octet matrix element 
(0g^( 3 Si)) because its dependence is cancelled in the ratio <j sd /<j nd . So, measuring this 
ratio R(ip) can determine the gluon fraction f g precisely, provided that the color-octet gluon 
fragmentation is the dominant mechanism for the J/ip production at large p?- 

One more thing that must be noted in the above calculations of the ratio R(ip) = 
usd/^nd is the approximation of neglecting the fragmentation function smearing. If the 
Pt distributions of the SD and ND J/ip production are much different, the smearing effects 
will influence the ration R(ip) and the extraction of f g from this ratio. To see these effects, 
we calculate the ratio R(ip) = ctsd/^nd for different values of p^(mm), which is shown 
in Table I. From this table, we can see that the ratio R(ip) is almost a constant (with a 
fluctuation less than 10%) as a function of pr(mm) for both the hard quark and hard gluon 
cases. This implies that the ratio R(ip) is not sensitive to the smearing, and we may neglect 
the smearing effects in the calculations of the ratio R(ip) and the extraction of the gluon 
fraction f g . 

In Table I, we also give the result for the soft gluon parametrization. The ratio R(ip) for 
the soft gluon is larger than the hard gluon by a factor < 2 at p^imin) > lOGeV. 

Experimentally, the non-diffractive background to the diffractive J/ip production must 
be dropped out to obtain useful information of the above calculations. Theoretically, the 
SD J/ip production events can be distinguished from those non-diffractive events by per- 
forming the Rapid Gap (RG) analysis. However, the acceptance of the RG will affect this 
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analysis. Here, we adopt the existing results of the background estimate obtained by the 
CDF diffractive Dijet experiment, where they give the non-diffractive background to the SD 
evens to be 20% ||10|| . By the same reason, we expect that the non-diffractive background 



to the SD J/ip production is about 20%. 

Finally, we discuss the theoretical uncertainty coming from the choice of the factor D 
of Eq.(0). The factor D represents the momentum fraction of the pomeron carried by the 
hard partons with the standard pomeron flux. In our calculations, we use the renormalized 
factor, which is about 1/9 at the Tevatron energy region. But the value cited here is not 
unique, because it may change with different choices of the parameters such as Mq and 
^ m ax in Eq.(0). The momentum fraction D can be measured in the diffractive processes at 
various collider faculties. At the Tevatron, the CDF have determined the fraction D to be 



0.18 ±0.04 fLOfl , which is well below the range 0.4 < D < 1.6 reported by ZEUS. If we adopt 
the CDF measurement, there must be difference in the above calculations of the ratio R, 
which is also shown in Fig.6. The shaded region in Fig. 6 represent the range of the ratio R 
calculated as a function of f g by using the fraction D = 0.18 ± 0.04. 

As discussed in previous studies [H)|,|2"0"|] , color-octet mechanism is crucially important to 



direct J/ip production (excluding the contributions from b decays and the Xc and ip' feed- 
downs) at large pr, here the color-octet mechanism is also crucially important to the direct 
J/ip production in the diffraction region. If only considering the color-singlet contributions 
(mainly coming from gluon fragmentation) the SD direct J / ip production rate will be smaller 
than the curves shown in Fig.3 by a factor of 50. This will make the measurement of the 
SD direct J/ip production very difficult at present luminosity at the Tevatron. That is to 
say, the SD direct J/ip production can also be regarded as another important test for the 
color-octet mechanism. 

As a final remark, we note that our proposal, the diffractive J/ip production, may also 
be used to extract the gluon fraction in photoproduction at the HERA. There are more 
diffractive events at the HERA than that at the Tevatron. So, more interesting results may 
be obtained. The work along this way is in progress |23| . 
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In conclusion, in this paper we show that the diffractive J/ip production at large px 
is sensitive to the gluon fraction of the pomeron f g . The measurement of this process at 
the Tevatron would provide a determination of f g . We have also discussed the uncertainties 
caused by the renormalised factor D and the gluon parameterizations of the pomeron. These 
uncertainties, however, can be reduced by combining other experimental measurements such 
as dijet diffractive production and W diffractive production. We believe that with the 
proposed SD J / ip production presented in this paper, we will get a better understanding for 
the property of the pomeron. And also, the SD direct J ftp production will provide another 
crucial test for the color-octet production mechanism. 
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TABLES 



P T {mm)(GeV) 


5.0 


6.0 


7.0 


8.0 


9.0 


10.0 


11.0 


12.0 


13.0 


R{ip) for hard quark (%) 


0.17 


0.16 


0.15 


0.14 


0.14 


0.14 


0.14 


0.14 


0.14 


R{ip) for hard gluon (%) 


0.98 


0.93 


0.90 


0.86 


0.84 


0.85 


0.86 


0.88 


0.87 



TABLE I. The ratio R(ip) = ctsd/o'nd as a function of pr(miri). 
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FIGURES 
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Fig.l Sketch diagram for the SD J ftp production by Pomeron exchange at 
the Tevatron. 
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Fig.3 The SD J/\|/ production cross section as a function 
of the minimum P T (\|/) for 2.0<r|(\|/)<4.0. 
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Fig.4 The diffrential cross section ove P T for the SD J/\|/ as a 
function of P T for the hard gluon (dashed line) and soft gluon 
(dotted line) parametrizations in the pomeron. 
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Fig.5 The ND forward J/\|/ production cross section 
for 2.0<r|(\|/)<4.0. 
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Fig.6. The ratio of R(\)/)=o SD /o ND vs the gluon fraction of the 

pomeron f . The solid line corresponds to the momentum fra- 

ction D used as the renormalized factor as in Eq.(7). The sha- 
ded region represent the range of the ratio R limited by the 
measured fraction D=0.18±0.04 by the CDF. 
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